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We study the nuclear modification of angular and momentum correlations between heavy quark 
pairs in ultra-relativistic heavy-ion collisions. The evolution of heavy quarks inside the thermalized 
medium is described via a modified Langevin approach that incorporates both elastic and inelastic 
interactions with the medium constituents. The spacetime evolution of the fireball is obtained 
from a (2-|-l)-dimensional viscous hydrodynamics simulation. The hadronization of heavy quarks 
is performed utilizing a hybrid model of fragmentation and coalescence. Our results show that the 
nuclear modification of the transverse momentum imbalance of DD pairs reflects the total energy 
loss experieirced by the heavy quarks and may help us probe specific regions of the medium. The 
angular correlation of heavy flavor pairs, especially in the low to intermediate transverse momentum 
regime, is sensitive to the detailed energy loss mechanism of heavy quarks inside the QGP. 


I. INTRODUCTION 


It has been shown that a new state of matter, com¬ 
monly referred to as the strongly interacting quark-gluon 
plasma (sQGP), can be produced in ultra-relativistic 
nucleus-nucleus collisions, such as those performed at the 
CERN Large Hadron Collider (LHC) and the BNL Rel¬ 
ativistic Heavy-ion Collider (RHIC). This highly-excited 
nuclear matter exhibits many remarkable properties [H- 
[^, among most important observations is the strong 
anisotropic collective flow exhibited by the hot and dense 
fireball and the final emitted hadrons , which can 

be nicely described by relativistic hydrodynamics sim¬ 
ulations [iM3- Another result is the strong modifica¬ 
tion of large transverse momentum (px) hadrons and jets 
as compared to the reference of proton-proton collisions 
[iMl- Many high-px phenomena and jet-related ob¬ 
servables can be well understood from parton energy loss 
which originates from the interaction between the propa¬ 
gating hard partons and the hot and dense QCD matter 

Heavy quarks and heavy flavor mesons are of partic¬ 
ular importance. Because of their large masses, heavy 
quarks are expected to lose less energy than light flavor 
partons due to the so-called “dead cone” effect, i.e., the 
phase space of collinear radiation is s upp ressed by the 
kinematics compared to light partons j22l |. In contrast, 
experimental measurements have shown that the nuclear 
modification factors Raa and the anisotropy coefficients 
V2 for heavy flavor mesons are comparable to light flavor 
hadrons [231 - 1^ . 

Many phenomenological models have been developed 
to investigate the propagation and energy loss of heavy 
^arks inside the hot and dense nuclear medium |29l - 
l3^ . In our earlier work we have developed a 

transport model to simulate the evolution dynamics of 
heavy flavors in relativistic heavy-ion collisions. Both 


collisional and radiative energy losses of heavy quarks 
are included in our model via an improved Langevin ap¬ 
proach, where the gluon radiation is treated as a recoil 
force exerted on the propagating heavy quarks. The re¬ 
coil force is obtained from the medium-induced gluon ra¬ 
diation spectrum as calculated from perturbative QCD 
energy loss formalisms. The hadronization process of 
heavy quarks after passing through the QGP is simu¬ 
lated utilizing a hybrid model of fragmentation and co¬ 
alescence, where the momentum dependence of the rela¬ 
tive probability of the two hadronization mechanisms is 
calculated according to the Wigner functions of the heavy 
hadrons. The interactions of heavy flavor mesons with 
the hadron gas were also incorporated in Ref. [s^ by 
utilizing the Ultra-relativistic Quantum Molecular Dy¬ 
namic Model (UrQMD) [i^. With the above setup, we 
can provide a reasonable description of the nuclear mod¬ 
ification of heavy flavor mesons at the LHC and RHIC 


In addition to single inclusive observables, correlations 
of final state particles have played a very important role 
in the study of jet quenching and parton energy loss. 
In the light flavor sector, there have been a wealth of 
experimental measurements and theoretical studies on fi¬ 
nal state correlations, such as back-to-back dihadron and 
photon-hadron correlations [dll - ld^ . as well as the mo¬ 
mentum imbalance of back-to-back jet pairs and photon- 
jet pairs (ttI - Is^ . In the heavy flavor sector, the angular 
correlations of heavy flavor pairs have been studied in 
Ref. [U, [ 534 ^ , and the momentum imbalance of heavy 
flavor pairs has been investigated in Ref. [fi^ [^. In 
this article, we present a detailed study of both angular 
and momentum correlations for heavy flavor pairs u sing 
our heavy flavor transport model built in Ref. [3511371 - 
l39j| . The dependence on the collisions centrality and the 
transverse momentum cuts are investigated for both an¬ 
gular and momentum correlations. We will show that 
while the momentum imbalance of back-to-back heavy 
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flavor meson pairs is mostly sensitive to the total en¬ 
ergy loss of heavy quarks in the dense QGP medium, 
the angular correlations probe the transverse momentum 
broadening of heavy quark propagation and are sensitive 
to different parton energy loss mechanisms, especially in 
the low transverse momentum region. 

This article is organized as follows. In Sec|lTl we 
present our theoretical framework for simulating the in¬ 
medium evolution and hadronization of heavy quarks in 
relativistic heavy-ion collisions. In Sec lIIIl we provide 
a detailed investigation on the angular correlations of 
heavy flavor pairs and show their sensitivity to different 
energy loss mechanisms. The transverse momentum im¬ 
balance between heavy flavor pairs is presented in SeclIVI 
Sec rVI contains our conclusions. 


II. HEAVY FLAVOR DYNAMICS 

In the limit of multiple scatterings, the momentum 
exchange of heavy quarks in each scattering with the 
medium constituents is small. In this case, the evolu¬ 
tion of heavy quarks inside a thermalized medium via 
multiple quasi-elastic scatterings can be described by the 
Langevin equation. However, as one moves to the high px 
region, quasi-elastic scatterings alone are no longer suf¬ 
ficient to describe the heavy quark evolution and energy 
loss in hot and dense nuclear medium, and the contri¬ 
bution from medium-induced gluon radiation has to be 
included. To incorporate both ener gy l oss processes, we 
follow our previous framework [ssl. l36l| and modify the 
classical Langevin equation as follows: 

^ = -VD{p)p + i + fg- ( 1 ) 

In the above equation, the first two terms on the right 
hand side are the drag force and the thermal random 
force experienced by heavy quarks when they scatter 
with the medium constituents; they are inherited from 
the classical Langevin equation. In this work, we assume 
the fluctuation-dissipation theorem still holds, rjjyiji) = 
k/{2TE), where k is the momentum space diffusion coef¬ 
ficient of heavy quark defined as — 

t'). On the right hand of the above equation, a third term 
fg is introduced to describe the recoil force experienced 
by heavy quarks when they radiate gluons. The proba¬ 
bility for a heavy quark to emit a gluon within the time 
interval [t,t + At] is evaluated according to the average 
number of radiated gluons: 

Prad(i, At) = (iVg(t, At)) = J dxdk]_ (2) 

where At is chosen to be sufficiently small to ensure 
< 1- dNg/dxdk'j^dt is the gluon radiation 
spectrum and in our work is obtained from the higher- 


twist energy loss calculation [63 - 1^ : 

dNg _ 2asP{x)q . 2 / 1 - t A / fcj 

dxdk^dt 7rfc^ ^ 2 t^ / \k^ + x'^M'^) 

(3) 

In the equation, x is the fractional energy carried by 
the radiated gluon, k± is its transverse momentum, 
Tf = 2Ex{l — x)/{k\_ -I- x^M^) is its formation time and 
P{x) is the vacuum splitting function, q is known as 
the gluon transport coefficient, which is related to n 
via q = 2 kCa/Cf- With this setup, there is only one 
free parameter in our model. In the literature, one of¬ 
ten quotes the spatial diffusion coefficient D of heavy 
quark to characterize its interaction strength with the 
medium; it is related to the momentum drag and dif¬ 
fusion as: D = r/[M77xi(0 )] = 2T^/k. We will follow 
such practice and quote the values of D in the following 
discussions. 

Using our transport model, we can simulate the propa¬ 
gation of heavy quarks inside the QGP matter created in 
relativistic heavy-ion collisions. The spacetime evolution 
of the local temperature and fluid velocity of the dense 
nuclear matter is simulated with a (2-|-l)-dimensional 
viscous hydrodyuamic model lu each time step, 

we boost the heavy quark into the local rest frame of 
the fluid cell in which the energy and momentum of a 
given heavy quark is updated according to our modified 
Langevin equation before it is boosted back to the global 
(laboratory) frame. 

Before the start of the hydrodynamical evolution (at 
To = 0.6 fm/c), the initial entropy density is obtained 
via a Monte-Garlo (MG) Glauber model. The same MG 
Glauber model is utilized to initialize the spatial dis¬ 
tribution of initial produced heavy quarks. The initial 
momentum space distribution of heavy quarks is calcu¬ 
lated in leading order perturbative QCD when calculat¬ 
ing the production and nuclear modification of single in¬ 
clusive heavy flavor mesons. The nuclear shadowing ef¬ 
fect in nucleus-nucleus collisions has been included by 
utilizing the EPS09 parametrization [i^. In the study of 
heavy flavor correlations, the initial angular and trans¬ 
verse momentum distributions of heavy flavor pairs are 
obtained from Pythia 8 simulation. The evolution of 
heavy quarks inside the QGP ends when they reach fluid 
cells with local temperature below (165 MeV). On 
the freeze-out hypersurface, heavy quarks hadronize into 
color neutral bound states using our hybrid model of frag - 
mentation and coalescence developed in Ref [H, ISq . 
In such a hybrid model, high momentum heavy quarks 
tend to fragment directly into hadrons, but low momen¬ 
tum heavy quarks are more probable to combine with 
thermal partons from the medium to form heavy flavor 
mesons. The relative probability between these two pro¬ 
cesses is determined by the Wigner functions based on an 
instantaneous coalescence model . Within our coales¬ 
cence model, the momentum distribution of the produced 
hadron can be calculated once a heavy quark is chosen 
to recombine with a thermal parton. If the heavy quark 
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FIG. 1: (Color online) D meson nuclear modification factor 
Raa in (a) central Pb-Pb collisions and (b) central Au-Au 
collisions. Different energy loss mechanisms are used for com¬ 
parison. 


is selected to fragment, Pythia 6 [t^ is used to simulate 
the fragmentation process using Peterson fragmentation 
function. The rescattering of the produced heavy mesons 
in a hadron gas will not be included in this work, since 
our previous study [s^ has shown that its effect can be 
reasonably approximated by slightly increasing the mo¬ 
mentum space transport coefficient of heavy quarks in 
the QGP phase. 

In our model, the transport coefficient of heavy quarks 
is determined by comparing to the experimental data on 
heavy meson nuclear modification at high px and then 
applied to calculations at lower px and other heavy flavor 
observables. In Fig. [U we provide calculations of the D 
meson Raa in 2.76 TeV central Pb-Pb collisions and 200 
GeV Au-Au collisions. Different energy loss mechanisms 
are utilized and compared to the experimental data. We 
can see that with the simultaneous inclusion of collisional 
and radiative processes, the use oi D = 5/(27rr) provides 
the best description for the px dependence of the D me¬ 
son Raa- If only quasi-elastic scattering or only gluon 
radiation are included, one can also obtain a reasonable 
description of the experimental data if the diffusion co- 
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FIG. 2: (Golor online) Angular correlations of cc pairs in 
central Au-Au collisions. Here the LO pQCD back-to-back 
approximation is used for cc pair initialization. 
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efficient is adjusted to D = 1.5/{2 ttT) ot D = 3/{2ttT), 
respectively. The uncertainty introduced by a different 
velocity dependence of the longitudinal and transverse 
transport coefficients (kl and kx) for the collisional en¬ 
ergy loss (denoted by the star notation) is investigated as 
well: the velocity dependence calculated to the leading 
logarithm is taken from Ref. [2^. In Fig. [T]we show that 
an adjustment of the diffusion constant to D = 2/(27rT) 
is able to reproduce our earlier results obtained with 
velocity-independent transport coefficients. Fig. [T] im¬ 
plies that the single particle spectrum appears not suffi¬ 
cient to distinguish between the contributions from dif¬ 
ferent energy loss mechanisms. In the next section, we 
will show that this large uncertainty can be reduced by 
investigating the correlations of heavy flavor pairs. In 
this work we perform the calculations for 200 GeV Au- 
Au collisions at RHIG. 


III. ANGULAR CORRELATIONS OF HEAVY 
FLAVOR PAIRS 

In this section, we study the angular correlations of 
heavy flavor pairs and their nuclear modifications in rel¬ 
ativistic heavy-ion collisions. The dependence of the nu¬ 
clear modification effect on the centrality cut and the 
transverse momentum px cuts will be investigated in 
some details. 

As a first academic study, we initialize heavy quark 
pair production using LO pQGD calculation, where c 
and c pairs are assumed to be back-to-back and have 
the same magnitudes of transvese momentum. We cal¬ 
culate the angular correlation functions of cc pairs after 
they traverse a realistic QGP medium created in central 
Au-Au collisions at RHIG. The result is shown in Fig. 
[21 in which Vp^ir denotes the number of cc pairs (later 
for D — D pairs within the selected trigger class) so that 
our correlation function is normalized to unity. The same 
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values of transport coefficients are used here as in Fig. [T] 
We can see that while different energy loss mechanisms 
provide similar amount of suppression of single D mesons, 
they result in very different angular correlation functions 
of cc pairs in the final state. Since the medium-induced 
gluon radiation is mostly dominated by the transverse 
momentum kicks exerted on the radiated gluons, the ef¬ 
fect on the propagation directions of the heavy quarks is 
very mild. One can see that the angular distribution of 
the final state cc pairs still peaks around = tt if only 
radiative energy loss is included. However, quasi-elastic 
scatterings are very effective in changing the heavy quark 
propagation directions, and the final angular distribution 
is widely spread out. Such result implies that the angular 
correlations of heavy flavor pairs provide a possibility to 
distinguish different energy loss mechanisms. 

One interesting observation in Fig. [2] is that purely 
collisional energy loss leads to a peak around Ai/) = 0 in 
the angular distribution. This indicates that a lot of low 
momentum cc pairs move collinearly in the final state 
even though they are initialized as back-to-back pairs. 
This is caused by the strong radial flow of the hydro- 
dynamic background which boosts the cc pairs into its 
own expanding direction. Such effect was historically re¬ 
ferred to as the “partonic wind effect” in Ref. [^; it 
might lead to some enhancement of quarkonium regener¬ 
ation in a QGP matter. We verify this effect by solving 
the Langevin equation in the global (laboratory) frame, 
i.e., the medium flow effect is shut off manually and only 
the temperature profile affects the heavy quark evolu¬ 
tion. With the flow effect turned off, the peak around 
= 0 disappears and the angular distribution of cc 
pairs becomes flat. The flat distribution also means that 
the final charm quark propagation directions can be en¬ 
tirely randomized by the strong quasielastic scatterings 
inside the QGP. In Fig. [2l we also observe that collisional 
energy loss with both velocity independent and velocity 
dependent kl and kt lead to similar angular correlation 
functions as long as the diffusion coefficients are properly 
adjusted to describe the observed D meson i?AA as shown 
in Fig. [TJ In the rest of this work, we will adopt the ve¬ 
locity independent k as assumed in our earlier study [3^ . 

The above study of the correlations between cc pairs is 
only for academic purposes. For the final state correla¬ 
tions of DD pairs, we first note that there might be more 
than one cc pair produced in one single event of heavy- 
ion collisions (which depends on the kinematic cuts and 
collision centrality). To take this effect into account, for 
each event we pair every D meson with all D mesons, and 
analyze all final state DD pairs. It is clear that there are 
uncorrelated DD pairs (background) contributing to the 
final correlations, especially with lowpx cuts and in more 
central collisions. 

For a more realistic simulation, Pythia 8 is utilized 
to simulate the initial cc pair production. We adopt the 
tunings of the Pythia 8 generator 0 that reproduce 
the spectrum of single charm hadrons in proton-proton 
collisions at RHIG as measured by the STAR Collabora- 



FIG. 3: (Color online) (a) px distribution of single inclusive 
D mesons, and (b) the angular correlations of DD pair in 
proton-proton collisions at RHIC obtained from Pythia 8 
simulation. 


tion [t^ [see Fig. 3(a)| , and then use the same parameter 
setting to obtain the angular correlation functions o f DD 
pairs in proton-proton collisions as shown in Fig. |3(b)[ 
We can see that the DD angular correlations are pretty 
flat if no transverse momentum cut is applied. If one im¬ 
poses a 2 GeV px cut for the trigger D or D meson, a peak 
appears around A(f> = tt. With higher and higher px cuts 
applied for the trigger and associated heavy mesons, the 
away side peak becomes more and more prominent. An¬ 
other interesting observation is that with a high enough 
px cut, a near-side peak begin to appear; this is the con¬ 
tribution from the initial gluon splitting process. The 
above obtained results serve as a baseline for the fol¬ 
lowing study of DD correlations in Au-Au collisions at 
RHIG. 

In Fig. 131 we show the angular correlation functions of 
DD pairs for 200 GeV central Au-Au collisions at RHIG. 
Fig. 4(a) shows the result when a 2 GeV px cut is ap¬ 
plied to the trigger D/D meson. Similar to cc pair cor¬ 
relation shown in Fig. [2l the purely radiative energy 
loss does not strongly affect the angular correlation very 
much [comparing the dashed line in Fig. |4(a)| and the 
solid line in Fig. |3(b)|. However, purely collisional en- 
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FIG. 4: (Color online) D-D angular correlations in central 
Au-Au collisions at RHIC. Different energy loss mechanisms 
and different pt cuts are shown for comparison. 


ergy loss is very effective in suppressing the away-side 
peak. Elastic scatterings also produce a near side peak 
around A(j) = 0, this effect diminishes when the coupling 
between the heavy quarks and the collective flow of the 
QGP is switched off. The realistic scenario in which a 
combination of both energy loss mechanisms is enabled 
falls between the above two curves. While these findings 
may help distinguish between different energy loss mech¬ 
anisms, note that systematic uncertainties exist for the 


radiative energy loss at low px since currently the de¬ 
tailed balance between gluon emission and absorption is 
not implemented and we employ a low momentum cut¬ 
off below which heavy quarks can only interact via elas¬ 
tic scattering to retain detailed balance for heavy quarks 
that have thermalized with the bulk QCD medium. An¬ 
other feature of DD correlations is the large background 
contributed from uncorrelated DD pairs, as compared to 
cc angular correlations. In our calculation, there are on 
average around 10 pairs of cc’s produced in each central 
Au-Au collision; the background contribution diminishes 
when one moves to peripheral collisions or larger px cuts 
are applied. For exam ple, (A^pair) drops below 1 per trig¬ 
ger event in Fig. 4(a) when the centrality is greater than 
about 50%. Fig. 4(b) shows the same results, but with 
the application of a 4 GeV px cut for the trigger D ox D 
meson and a 2 GeV cut for the associated meson. One 
can see that the difference between collisional and radia¬ 


tive energy losses becomes smaller than that in Fig. 4(a) 


The difference almost disappears with higher enough px 
cuts applied, see e.g.. Fig. 4(c)[ As shown by the initial 
spectrum of charm quarks at RHIG energy [Fig. 3(a)|, 
with the Px cuts applied in Fig. |4(b)| and 4(c)[ on av¬ 
erage there should be less than one pair of cc produced 
in each event, and thus very few uncorrelated DD pairs 
contribute to these two scenarios. 


To further quantify the medium modification of the 
angular correlations of heavy flavor meson pairs, we cal¬ 
culate the variances of the angular distributions for both 
near side (— 7r/2,7r/2) and away side (7r/2, 37r/2). The 
result is shown in Fig. [5] where the variances are plot¬ 
ted as functions of the participant numbers fVpart in Au- 
Au collisions at RHIG. The points at Apart = 0 cor¬ 
respond to the baseline proton-proton collisions. Note 
that the variance of A(j) is defined as {{A(j) — (Ai^))^) = 
{A(j)^) — {A(j))'^. And for the purpose of a better il¬ 
lustration, we have drawn a horizontal dashed line at 
7r/-\/l2 which is the variance for a uniform distribution 
in (— 7r/2,7r/2) or in (7r/2,37r/2); a smaller variance than 
7r/-\/l2 implies a peak around Acj) = 0 or = tt, and a 
larger variance than t: l\/V2 implies a dip on the near or 
away side. 

Figure [5 (a) | shows the variances for the near side cor¬ 
relations. When a 2 GeV cut is applied to the trigger 
Px (the upper panel), there is a dip structure in the 
DD correlations in proton-proton collisions and periph¬ 
eral Au-Au collisions (therefore, the variance is larger 
than the value of a uniform distribution). With purely 
radiative energy loss, such dip gradually approaches a 
uniform distribution as one moves to more central colli¬ 
sions. However, with purely collisional energy loss, the 
variance first drops below 7r/-\/T2 when Apart is not very 
large, and then increases and gradually approaches to 
the uniform distribution value as one moves to more cen¬ 
tral collisions. The smaller variance than 7r/-\/T2 in cen¬ 
tral and mid-central collisions indicates there is a peak 
formed around A(j) = 0 due to the “partonic wind effect” 
as discussed earlier. The non-monotonic behavior of the 
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FIG. 5: (Color online) Variances of (a) near-side and (b) 
away-side distribution in DD angular correlation function as 
functions of the participant number. Different energy loss 
mechanisms and different px cuts are shown for comparison. 


starts from a value below 7r/-\/T2 for proton-proton colli¬ 
sions and decreases with the increase of A^part in An-An 
collisions. This behavior can be understood from the in¬ 
medium energy loss of heavy quarks as follows. The final 
D/D mesons within certain px cuts originate from ini¬ 
tial c/c quark which have much higher momentum and 
thus smaller variances for their angular distributions [see 
Fig. |3(b)| . Since the angular correlations do not change 
very much at highpx after in-medium evolution [see e.g., 
Fig. |4(c)| , one obtains more peaked angular distribution 
for heavy flavor pairs for more central collisions. We 
further verify this by applying the transverse momen¬ 
tum cuts to the initial heavy quarks instead of the fi¬ 
nal observed heavy mesons. For such fixed sample of 
heavy quark/meson pairs, we do find that the variance 
of the angular distribution increases monotonically to¬ 
wards TT/\/\2 when increasing Apart- Again, we find that 
the difference among different energy loss mechanisms is 
larger for low px meson pairs and becomes smaller for 
larger px meson pairs. 

The variances for the away-side angular distribution 
are shown in Fig. |5(b)| With small px cut (2 GeV), 
purely radiative energy loss smoothes the away-side peak 
in DD correlation distribution and the variance mono¬ 
tonically approaches towards the uniform distribution 
value. However, with purely collisional energy loss, the 
variance first increases above the uniform distribution 
limit which indicates a depletion of DD pairs around 
A<^ = TT, and then decreases towards Tr/y/^ due to the 
smoothing of the background contribution. With higher 
Px cuts, there is the interplay between two combina¬ 
tional effects: the transverse momentum broadening and 
the energy loss of heavy quarks inside the QGP. The 
former increases the variance towards the uniform dis¬ 
tribution limit and is more efficient for relatively lower 
Px, while the latter enhances the contribution from the 
charm quarks with higher initial px which reduces the 
variance of the final angular distribution. As a conse¬ 
quence, a non-monotonic behavior of the variance is ob¬ 
served as a function of collision centrality. In the next 
section, we will disentangle these two contributions - mo¬ 
mentum broadening and energy loss ~ to the value of the 
variance by utilizing the momentum imbalance of DD 
pairs. An alternative way to distinguish these two con¬ 
tributions - using the half width at the half maximum ~ 
was discussed in Ref. [s^. 

IV. MOMENTUM IMBALANCE OF DD PAIRS 


variance as a function of Apart originates from the in¬ 
creasing contribution from the uncorrelated DD pairs in 
more central collisions which smoothes out the angular 
distribution and make the variance approach the uniform 
distribution value. 


With the application of larger px cuts [the middle and 
lower panels in Fig. 5(a)|, the variance at the near side 


In addition to the nuclear modification of heavy ffavor 
pair angular correlations, another interesting observable 
for studying the correlations of heavy mesons is their 
momentum imbalance. For this study, we may define 
the transverse momentum ratio: = px,asso/PT,trig for 

DD pairs, where Px.trig and px,asso are the transverse 
momenta of the leading and subleading mesons, whose 
azimuthal angles are denoted as ^trig and (j)a.sso- For each 
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FIG. 6: (Color online) DD transverse momentum imbalance 
distribution for different centralities and different px cuts. 


event, we select the D or D meson that has the high¬ 
est transverse momentum as the leading (trigger) meson. 
On the back side, we look for the anti-particle {D or D 
meson) with the highest transverse momentum and se¬ 
lect it as the subleading (associated) meson. An angular 
cut is also applied for the away-side subleading meson: 

I<('asso — (/i’trig I > 27r/3; this should be helpful to reduce the 
background contribution and select more truly correlated 
pairs. 

In Fig. |6l we show the event distribution (normalized 
to one D 01 D trigger) of the px ratio (xt) of DD pairs for 
different centrality bins and kinematic cuts. The results 
are shown for two different px cuts: px,trig > 4 GeV, 
PT.asso > 2 GeV in Fig. Ela), and Px.trig > 8 GeV, 
PT.asso > 4 GeV in Fig. (IDb). Two observations can be 
found in Fig. [51 as one moves from proton-proton col¬ 
lisions to peripheral and subsequently to central Au-Au 
collisions, (1) a smaller number of DD pairs per triggered 
event can be seen, and (2) the xx distributions for both 
sets of px selections tend to shift to the region of smaller 
xx (i.e., larger momentum imbalance). These are both 
due to stronger energy loss of heavy quarks inside the 
dense nuclear medium in more central collisions. As dis¬ 
cussed earlier, there should be very few uncorrelated DD 



part 

FIG. 7: (Golor online) The average values of DD momentum 
imbalance distribution as a function of the participant number 
for different px cuts. 


pairs contributing to the correlation functions here with 
the chosen px cuts. In the subfigures of Fig. [51 we also 
present the ratios between nucleus-nucleus collisions and 
proton-proton collisions, i.e., Iaa as defined by 


1 dNpa^ij^ 


1 

1 

Ntiig 

dx \pp 


One can see that the values of Iaa are above unity at 
small xx and below at large xx. At large xx, it decreases 
with the centrality of the collisions. 

To have a cleaner description of the transverse mo¬ 
mentum imbalance of DD pairs as a function of collision 
centrality, we show in Fig. [7] the average value of the 
Px ratio (or xx) between the subleading and the leading 
mesons as a function of Vpart for the two sets of kine¬ 
matic cuts, where Vpart = 0 corresponds to the proton- 
proton reference. One can see that as one moves from 
proton-proton collisions to central Au-Au collisions, the 
transverse momentum imbalance of DD pairs increases 
(xx value changes about 6%). It might be easier to ob¬ 
serve a larger medium modification to the momentum 
imbalance of DD pairs at the LHG with even higher px 
cuts, as shown in Ref. [g^. 

With the knowledge of xx one may further classify 
the events within the same centrality bin. It has been 
shown in Ref. [52| that for y-jet larger xx corresponds to 
events in which hard partons are produced at the edges 
of the QGP fireballs while smaller xx values are con¬ 
tributed by those pairs produced at the centers. In Fig. 
[8] we implement a similar study to investigate the corre¬ 
lation between the xx values of DD pairs and the pro¬ 
duction positions of the initial cc pairs. We observe that 
for heavy flavor pairs, smaller xx corresponds to events 
initially produced at the edge of the QGP fireballs in 
which one heavy quark travels outside the medium with¬ 
out much interaction while its partner traverse the whole 
QGP hreball and loses significantly amount of energy. 
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FIG. 8: (Color online) The density distribution of the initial 
cc production points (sinit, 2/init) when the triggered D or D’s 
are taken along the out-of-plane directions (|<(>trig — 7r/2| < 
7r/6) for 0-10% Au-Au collisions at 200 GeV (pT.trig > 4 GeV 
and pT,asso > 2 GeV). Different values of xt of the final DD 
pairs are compared: (a) xt € [0.2, 0.4], (b) xt € [0.4, 0.6], (c) 
Xt G [0.6, 0.8], and (d) xt € [0.8, 1.0]. 


For larger xt values, initial cc pairs are more likely to 
spread out smoothly over the QGP. Our study thus con¬ 
firms the pattern found in (Fig. 9). This allows us 
to utilize this momentum imbalance to probe different 
regions of the dense nuclear matter and gain knowledge 
on the path length dependence of parton energy loss as 
well. Note that in Fig. IHIHI we present our calculations 
with both collisional and radiative energy loss switched 
on. We have verified that since the momentum imbalance 
is directly related to the total amount of energy loss of 
heavy quarks, its value does not strongly depend on the 
detailed energy loss mechanism as long as the transport 
coefficient is properly adjusted to describe the D meson 
Raa- 

The momentum imbalance allows us to refine the pre¬ 
vious study on the angular correlation function in Sec. 
uni In order investigate the competition between the two 
ingredients that affect the variance of the angular corre¬ 
lation function - momentum broadening and energy loss 
- as shown in Fig. [Sj the values of the variance are now 
calculated for different xt bins. Here we utilize the away- 
side variance obtained with the intermediate pT triggers 
as an example, which has the best potential to constrain 
the different energy loss models. In Fig. [Hlwe can observe 
when Xt is small (< 0.4), it is easier for the large energy 
loss to overwhelm the effect of the momentum broaden¬ 
ing. With the pure radiative energy loss, the away side 
variance decreases monotonically; to the contrary, the 



part 

FIG. 9: (Color online) The xt dependence of the away-side 
variances of the mid-px DD angular correlation functions. 
Different energy loss mechanisms are shown for comparison. 

pure collisional energy loss would increase this variance 
in more central collisions since it smears the direction of 
heavy quarks more effectively inside the QGP; and the 
combined energy loss mechanism leads to a relatively flat 
variance when the participant number is large. However, 
in higher xt bins (or with smaller energy loss), the ef¬ 
fect of momentum broadening starts to dominate when 
Npa,rt is large and thus increases the away side variance 
towards the value of a uniform distribution. And the 
turning point between the decrease and increase of the 
variance starts earlier when a higher the xt cut is applied. 
Thus, the combination of the momentum imbalance and 
the angular correlation function may provide us a quan¬ 
titative tool to study the detailed energy loss dynamics 
of heavy quarks inside the QGP. 


V. CONCLUSIONS 

In this work, we have studied the correlations of heavy 
meson pairs in relativistic heavy-ion collisions. The 
transport of open heavy quarks inside a dense nuclear 
matter is simulated utilizing our modified Langevin equa¬ 
tion that incorporates both quasi-elastic scattering and 
medium-induced gluon radiation. The spacetime evo¬ 
lution of QGP fireball produced in Au-Au collisions at 
RHIG is simulated with a (2-|-l)-dimensional viscous hy¬ 
drodynamic model. A hybrid model of fragmentation and 
heavy-light quark coalescence is applied to simulate the 
hadronization of heavy quarks into their mesonic states. 

We have explored two relatively new observables for 




































9 


the nuclear modification of heavy flavors: the angular 
correlations and the transverse momentum imbalance of 
heavy meson pairs. Our results have shown that the 
transverse momentum imbalance of heavy meson pairs 
provides another measure of the energy loss of heavy 
quarks inside the QGP. Compared to the nuclear mod¬ 
ification of single inclusive meson spectrum, one of the 
advantages of this observable is that it is not only able 
to quantify the total energy loss of the heavy quark, but 
also provides an opportunity to help probe different areas 
of the QGP fireball. It is found that with the increase of 
collision centrality, the average px ratio between the sub¬ 
leading and the leading heavy flavor mesons decreases, 
which indicates a larger imbalance caused by the larger 
energy loss of heavy quarks. One the other hand, the 
angular correlation of DD pairs provides a good candi¬ 
date to characterize the momentum broadening of heavy 
quarks inside the QGP, and is also sensitive to differ¬ 
ent energy loss mechanisms. We find that while purely 
radiative energy loss does not affect very much the angu¬ 
lar distribution between DD pairs, the collisional energy 
loss is very effective in randomizing the propagation di¬ 
rection of heavy quarks, especially at low to intermediate 
Pt- To better quantify the medium modification of the 
heavy flavor pair correlations, we have also calculated 
the variances for both near-side and away-side DD an¬ 
gular distributions and investigate them as a function of 
collision centrality and transverse momentum cuts. By 
utilizing different bins of the momentum imbalance, the 
competition between the two effects - momentum broad¬ 
ening and energy loss - on the medium modification of 
the variance can be clearly seen on the away side when 


intermediate px cuts are implemented. 

Our study constitutes an important contribution to¬ 
wards a more quantitative understanding of the energy 
loss and momentum broadening of heavy quarks inside a 
hot and dense nuclear matter. If heavy flavor meson pair 
correlations are measured in the future, one can not only 
infer how much energy is lost from heavy quarks, but 
also gain some insights on the detailed interaction mech¬ 
anism between heavy quarks and the QGD medium; the 
later seems quite difficult by only studying the nuclear 
modification of single inclusive heavy flavor meson spec¬ 
tra. Nevertheless, it is worth noticing that in order to 
draw further quantitative conclusion on the heavy flavor 
dynamics from these correlation functions, it is crucial 
to develop a sound method for background subtraction, 
especially in the low px regime and in central collisions. 
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